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ABSTRACT 





These investigations generalize Miyata's synthesis of passive 
driving-point impedance functions by developing a step-by-step com- 
putational technique for augmenting a general driving-point impedance 
to insure that the real part of the augmented impedance is positive 


term by term. This goal has been accomplished and programmed under 





the provision that the real part of the original impedance is minimum 
reactive, i.e., has no zeros on the jw-axis. The latter requirement 


is not restrictive since zeros on the jw-axis can be removed a priori. 
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I. LNTRODUCTION 


>" method of transformerless synthesis provides a method 


The eae 
for separating a given impedance into additive components which indi- 
vidually satisfy certain special conditions that permit their realization 
in simple terms. This method of separation is based upon the even part 
of the given impedance and can be stated briefly as follows: Given a 
minimum reactive driving point impedance Z(s), find R(w) the real part 
of Z(s). Then split R(@) into a sum of even functions a) 2 (p20; 13..c ay 
each satisfying the requirements of Bee Ie Find AS) corresponding to 
each B®). The realization of Z(s) reduces to the realization of each of 
the a) and the series connection of the resulting one ports completes 
the synthesis. 

An inherent difficulty in the above technique lies in the require- 
ment that the R@) be positive for all w. To overcome this difficulty 
Miyata suggested the use of surplus factors. An example will illustrate 


this method. 


Let Z(s) be given as 


2 
Z(s) = $5. 70st (1.1) 
5 s +3s+4 
from which we obtain 
4 2 
R(w) = ee = R, @) + R.@) + R,.@) Cle2) 
Bee) l 0 


This so called "single-n split'’ method cannot be applied because 
R, ©) is negative. Therefore use a surplus factor (stl) to obtain an 


augmented impedance, Z_(s)- 


ie) 2 
Z_(s) = aa Z(s) = sie (1.3) 
s +4 s +7st4 





for which the numerator of ReZ_ (jw) contains only positive terms. 


6 2 
RGD bw _+30n +36 Cia 


2 
B' Ww) 

Obviously the use of surplus factors complicates matters since 
higher-order factors are often required. 


> 


At this point most authors considering the technique, including 
Miyata, generally abandon the use of surplus factors without discussing 
any general method for their selection and proceed to describe various 
groupings of the individual positive and negative Be which remain 
positive for all w. These grouping methods generally require factori- 
zation of the numerator with its attendant difficulties and no general 
method is available in the literature nor is it clear that such groupings 
are always possible. 

The purpose of this thesis is to investigate in detail the use of 
surplus factors to develop augmented impedances that are positive term 
by term. 

Chapter two is a general consideration of the Miyata method which 
describes the reasoning behind the technique and the mechanics of 
developing the electric circuit and component values from the original 
driving-point impedance when all terms in the numerator of ReZ(jw) are 


initially positive. The significant results of this section include: 


The circuit realization is transformerless, an important advantage 
over the Bott-Duffin technique. 


The complete circuit development consists of a simple unbalanced 
ladder network with series L's and shunt C's terminated in a resistance. 


By term grouping without factorization, the number of circuit elements 
required is less than for the modified Bott-Duffin for all driving-point 


(dp) impedance functions of second order or greater. The saving in 





number of elements becomes quite dramatic for higher-order polynomials 
so that even though the order of the original dp impedance must often 
be increased by the use of surplus factors, large savings are available. 

In general, chapter two is included to provide sufficient incentive 
to undertake the task of developing a general method of selecting 
surplus factors of minimum order. 

In chapter three, starting with a dp impedance that is the ratio of 
two general polynomials of the same order in s and a general augmenting 
polynomial (surplus factor), the real part of Z(s) is developed. Since 
the resultant coefficients of uae must be non-negative, these coefficients 
are the required constraints that must be satisfied to insure that ReZ(jw) 
is negative term by term. These coefficient constraints constitute a 
set of non-linear inequalities in the variables 6. (i=0,1,...,n) which 
are the coefficients of the augmenting polynomial (i.e., P(s)=6 +5. st...46 5°). 
Also in this section, the non-linear inequalities developed are linearized 
using a set of Sonendent transformations and solved by the Piylere method. 
The difficulties involved in the inverse transformation are circumvented by 
considering the form of the solved constraints which are the coefficients 
of the real part of the augmented dp impedance. The process is completed 
by demonstrating how a knowledge of the form of the solved constraints 
allows one to find the numerator of the real part by simple long division. 
Determination of the denominator is straightforward and the problem is 
solved. 

A digital computer program for part of the routine described above is 
developed in chapter four, and a sample problem worked in section five. 

In summary, a method has been developed which will invariable find 


an optimal augmenting polynomial such that the resultant augmented dp 





impedance has a real part which is positive term by term. The method 
is optimal in the sense that as many coefficients as possible of the 
real part are zero, thereby reducing the required number of circuit 
elements to a minimum. 

The value of this method lies in the fact that it will invariably 
produce a positive real part for the augmented impedance in a step-by- 
step fashion. In fact, the entire process could be programmed such that 
for any input minimum-reactive dp impedance the output would be the 
required circuit with accompanying element values. 

The single major disadvantage is that the augmented dp impedance may 
have a large number of terms which would consequently require a large 
number of circuit elements to realize (the circumstances that determine 


this possibility are discussed in chapter two). 
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II. GENERAL CONSIDERATION OF THE MIYATA SYNTHESIS TECHNIQUE 


The Miyata method for realization of a passive, minimum-reactive, 


> 


driving-point impedance proceeds as follows: For an impedance having 


no jw-axis zeros or poles written 


(2.1) 


where m denotes an even polynomial in s and n denotes an odd polynomial 
ii S.. 


Its real part for s=jw is 








n,m ane he) 6 
ReZ(jw) = = 5 Z (2720 
‘i 
or 
A +A 0+. +A A Aw Awe 
ae) 5 ——. Ow. 2 ne y) (2.3) 
BtB yw +...+B w Bw ) Bw ) Bw ) 


Each of the terms on the right-hand side, when the coefficients 
A Ayo e OAL are positive, may be shown to generate an impedance that 
can be realized without mutual inductive coupling. The series con- 
nections of such separate realizations yields the desired synthesis of 
Z(s). 

To see why this is so, consider the component impedances of the real 
part. The first term has n double-order zeros at s= infinity; the second 
has one double-order zero at s=0 and n-l zeros at s= infinity; the third 
term has two double-order zeros at s=0 and n-2 zeros at s= infinity; 
finally, the last term has all its n double-order zeros at s=0. There- 
fore, all of these terms have their double-order zeros located exclusively 
at s=0 and s= infinity. 
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When the real part of an impedance is zero at some jw-axis point, 
it need not necessarily follow that the entire impedance is zero. 
However, the jw-axis points under consideration are restricted to s=0 
and s= infinity and, since a minimum-reactive function must be purely 
real at these points, it follows that when the real part is zero the 
whole impedance must necessarily be zero. 

Next we use the namiiseane theory to show that these impedances are 
completely developable in simple terms. The Darlington theory shows that 
in a cascade development of an impedance each component lossless network 
produces one set of zeros out of the total complement of zeros which the 
even part of that impedance contains. Therefore, the network generated 
by the first cycle provides one set of zeros and leaves a remainder whose 
even part contains the remaining even-part zeros. The next network in 
the chain provides another set and leaves a second remainder with the 
leftover even-part zeros - and so forth. 

Applying this reasoning to the real part of the first term of the 
expansion of the real past of Z(s), we can predict that the corres- 
ponding impedances will have a zero at s= infinity and its reciprocal 
a pole at this point. Removal of this pole as a shunt capacitor is one 
Darlington cycle. The admittance remainder after this pole at s= infinity 
is removed has an even part that still has n-1 double-order zeros at 
s= infinity since so far only one has been used up. This admittance 
remainder must have a zero at s= infinity and its reciprocal must have 
a pole there. Removal of this pole as a series inductor constitutes the 
second cycle. The remainder now has an even part with n-2 double-order 
zeros at s= infinity and is zero at s= infinity. Its reciprocal must 


have a pole at this point and its removal as a shunt capacitance completes 
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the third cycle. This process continues yielding series inductive and 
shunt capacitive branches of total number n, the number of double-order 
zeros poSsessed by the original real part. When the last of these zeros 
has been accounted for, the remainder is reduced to a constant because 
its real part has no more zeros. The complete development consists of a 
simple unbalanced ladder network with series inductances and shunt cap- 
acitances terminated in a resistance. Subsequent terms are similarly 
developed. 

Using this single-n split technique, the total number of reactive 
elements is sie where n equals the highest degree of either the numerator 
or denominator polynomial in s. Comparison with the number of elements 
required for a Bott-Duffin synthesis is displayed in Table I along with 
the reductions due to the 1/2-n split technique which is subsequently 
discussed. The geometric character of the Bott-Duffin method compared 
with the arithmetic character of the Miyata method is strongly apparent 
in this table. Even though special devices are required to salvage the 
Miyata method when all the terms of the real part numerator are not 
initially positive, there appears to be considerable advantage to 
exploring the Miyata procedure where more elaborate impedance functions 
are involved. Before one gets too enthusiastic about the Miyata method, 
it would be well to consider the requirement that all the coefficients 
in the polynomial formed by the numerator of the real part of the 
driving-point impedance must be positive. 

A sufficient condition that the coefficients be positive can be 
obtained by recalling that a Hurwitz polynomial has positive coefficients. 
Thus if 


i ana +... tAQW (2.4) 
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has no zeros in the PRE ha fw plane, then A “9A are positive. 


or: 
Or if 


2 
A(-s ) = mym, - njn, (2:%5)) 


has no zeros in the left-half ag plane, then the A oA are positive. 


0°: 
Since the left-half se plane maps into a region of the s plane con- 
taining all points not more than plus or or 45° from the jw-axis, we 
can say the coefficients in eq. 2.4 are positive so long as there are no 
even-part zeros of Z(s) closer than 45° to the jw-axis. 

Although the coefficients in 2.4 may still be positive when some 
even-part zeros lie closer to the jw-axis than this sufficiency con- 
dition permits, the chances of this desired condition being fulfilled 
when even-part zeros are relatively close to the jw-axis become remote 
and vanish completely when a single pair of even-part zeros lie upon the 
jw-axis. 

In a manner of speaking, we can say that a polynomial that has some 
zeros above the 45° line in the s plane may still yield all positive 
coefficients if there are a sufficient number of zeros below the 45° 
lines to produce a compensating effect. Therefore, in a situation in 
which the even part of Z(s) has some negative coefficients, we can augment 
the real part with appropriately chosen factors to bring about the desired 
compensation. 

The method of compensation fails if the even part of Z(s) has any 
jw-axis zeros, since the number of necessary compensating factors tends 
toward infinity as one or more pairs of even-part zeros approach the 
jw-axis. With such a compensating scheme the number of terms in the real 


part and hence the number of necessary elements increases rapidly. 


14 





Orie method to reduce the number of elements is considered now, the 
1/2-n sige, °° Miyata points out that the detailed procedure for 
carrying out the basic method of decomposition may be varied by separating 
the real-part numerator polynomial into groups of terms rather than 


Separate terms. 


Consider 
(AHA ot. tA OK) CI Conn eT UTR 
ReZ( jm) = a ae 5 (2.6) 
B@ ) B(w ) 


where k=n/2 (mn even) k=(n+1)/2 (mn odd). 
>. i. nue who refers to this decomposition as the 1/2-n split, finds 
that it yields a smaller total number of elements than any other division 
of the numerator polynomial into parts. It requires fewer elements than 
the Bott-Duffin method for any degree n equal to or greater than 2. The 
impedance generated by the first term is developable into a simple unbalance 
ladder with series L's and shunt C's (total number k) terminated in a 
remainder whose even-part zeros are the same finite non-zero ones contained 
in the even part of this impedance. The second term is developable into an 
unbalanced ladder with series C's and shunt L's (total number n-ktl) and a 
terminal impedance whose even-part zeros are those that are contained in 
the even part of the second impedance. To the terminal impedances which 
appear as remainder functions we can apply a similar treatment and continue 
until the ultimate remainder functions are resistances. 

The process of augmentation by surplus factors is also necessary in 


order that Miyata's elegant method of determining an impedance from its 


real part may be used. The method proceeds as aenloses 
m, +n 
it I 
DD = a7) 
Mo ai no 
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and 


eae = Tee 


eZ 1 2 N(s 
EvZ(s) = oe (2.8) 
Ms = 2) 2 


Next consider an auxiliary function 


Z'(s) _ aii ae (2.9) 
mM, + ny ; ; 
Its even part is 
_ 
EvZ'(s) = sh ie “1 na (2.10 
m 2 

Set the two even parts to be equal; 

EvZ'(s) = EvZ(s) (2.1) 
which can be done provided 

m,'m, -- ny, My = 1, (2a) 
where the degree of m, ‘tn,’ is that of m,tn,- Polynomials my and 

ny can always be found to satisfy this requirements. Z'(s) however will 


have a numerator of degree two or more higher than its denominator and is 
therefore not positive real (p.r.). To find the p.r. function Z(s), carry 
out a long division on Z'(s) until the degree of the numerator does not 
exceed the degree of the denominator. That is, 


rot n,') m +n 


; = ip “ r r 
Z'(s) mn, +a, i q(s) +a Fa, rc 7 (2.13) 


i) 


where q(s) is the quotient polynomial and m tn is the remainder. 


Let zZ_(s) be 
mi + n 
z_(s) = oo (2.14) 
so that 
z(s) = Zs) + a(s) (2.15) 
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Recalling that EvZ'(s) = EvZ(s), 
EvZ(s) = EvZ (s) + Eva (s). (2.16) 
If the polynomial q(s) contains an even power of s then the right 
hand side of eq. 2.16 tends to infinity as s=jw tends to infinity. On 
the other hand, EV Eig) = R(w) is bounded at infinity. Therefore q(s) 
must be an odd polynomial and as a result 
EvZ _(s) = EvZ(s). (2717) 
Since the numerator of Z_(s) does not exceed the denominator Z_(s) 
has no pole at infinity or at any other finite frequency on the jw-axis 
since its denominator is that of Z(s). 
Briefly Zs) has no poles in the right-half plane or on the jw-axis, 
and its even part for s=jw is always non-negative. Therefore Zs) is 


positive real. To illustrate: 


Given ReZ (jw) -—— (2. ia) 


l1+wo 


2 
Z (s) = efi otis tt (2.17b) 
s +2s +2s +1 


from any method whatsoever. Suppose we want an impedance with the 
real part ReZ , (jw) = i iN which is 2.17a multiplied pomane and implies 
m=-s 

Then F(s) equals eq. 2.17b multiplied by a6 and the long division 


takes the form -2/3s 









4. 3 2 
3 2 -2/3s, - 4/3s. - s 


wee- + 2/3s 


and 


2 
z4(s) _ wos toahie 
s +2s + 2s + 1 
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Using this procedure any Z_(s) a ical fa) fie 2s. Cane be 


found from Zi (s). If we want an impedance with the real part 
Pa: 2 4 6 
ReZ(jw) = A, + A jw + Ajw + AW 


Z(s) may be generated as 
AZ. (s) + AZ, (s) + AyZ,(s) + A,2, (s) 
where continuing as above 


fi 2/3s° + 1/3s 
Z,(s) 3 2 
s + 2s + 2s +1 


. aa + nfope + 2/3s 
sae ie 
s +2s + 2s + l 


| 


Z,(s) is found by multiplying Z,(s) by wae and Z..(s) found by 


multiplying Z,(s) by e and performing the indicated long division. 
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III. GENERAL STEP-BY-STEP PROCEDURE FOR THE SELECTION OF SURPLUS FACTORS 


The topic of this chapter is the step-by-step procedure for deve lop- 
ing the real part of an augmented impedance that is positive term by term 
th : ae: ; al , 
from a general n -order, passive, minimum-reactive, driving-point 


impedance augmented by a general fugeerder augmenting polynomial. 


A. CONSTRAINT DEVELOPMENT 

This section is devoted to the development of the required constraints 
to insure a positive real part for the dp impedance and begins with a 
general impedance function and augmenting polynomial as follows: 


Consider first the general driving point impedance Z(s) 


N 
% as” N 
_ q-0 q_ _ N(s) 
Z(s) = N , 4 D(s) (el) 
Zed 
d= 


T 
P(is) = E 6 8 (3.2) 


The product of Z(s) and P(s)/P(s) is the augmented driving point 


impedance Z_ As); 


mM = 
os 
Q 
WY 
Q 


(3 .3) 


N 
ab) 
o~ 
¢) 
w 
Il 
eo a) 
2 2 
ry we 
tl aes i aes 
= 


NM =z 
“~~ 
co) 
9) 
wa) 
Se 


2 
IMA =@ 


me ao 


Se (3.4) 


May zl il 


ry Oo 





The even part of the augmented driving point impedance is 


m.‘'m -n.'n 
> 2 m2 NEMS 
EvZ_(s) =. 2 ; 2 = D'(s) (35) 
a oD | 
where 
(N+T) /2 (MET) /2 
T T 
2u 20 
a f= 
oT x Qo or ae? “ cee r° 
u=0 u=0 
r=0 r=0 
(3.6) 
(N+T-1)/2 (N+T-1)/2 
T T 
2utl 2utl 
ne i 
a5) lene =D 2s posers aa 
u= u=0 
r= = r=0 


and noninteger values of (NT)/2 and (N+T-1)/2 are rounded to the 
next lowest integer. 


Expanding the even part of Z_,(s) 


(M-T ) /2 (MT )/2 (N+T-1)/2 (NT-1)/2 
T T T T 
- u 2u Zur Zur 
Ev2(s) = = Mu eon Boner r® ~ ler z let r° 
u=0 =0 =0 =0 
r=0 r=0 r= r=0 
D'(s) (Sag) 
Combining products, 
(MT ) /2 (NFT-1)/2 (3.8) 
T Ab 
2 (utv ) 2(utv+1) 
EvZ (s) = ZX s Shbea> 2a) Ss § 6 
af ) re=0 ra ales ea EO vt l-r°2ubl-q rq 
q,r=0 r=0 
D'(s) 
Combining the difference, 
(MET) /2 
1 
2 (utv) 
EvZ9¢s) = > 6. as (of 3.29 
a (s) a (doy x? ung jv a en q? i 
r,q= Oa? 
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Since ReZ_ (jw) = Be) acu and the Miyata method requires that 
the real part be positive term by term, constraints must be developed 
to insure this requirement. First, consider the denominator D'(s) above 


D'(s) =m," - ny! = D(jw)D(-jo) = [D(ja) |” (3.10) 


which is never negative and is therefore of no concern. 


Thus the constraint is reduced to 
2 > << 
N'@ ) mm, ny 2, = 0 for O=»p=0 (3.11) 


The coefficient of each term in eq. 3.9 must be non-negative which 


requires that 


(N+T-1)/2 
(N+T) /2 
Gort2 


< 
ORO lege D eo ou nG ; Ogee led uae” uel ¢° — (2 xy 


and 


(N+T-2)/2 
(N+T)/2 
iE 
oF 
n=0 
u=0 
r ,q=0 


4mt&4 > 
: ra er ner j CH ieQu-r°2ut1-q" 5S (3.13) 


In equation 3.9 when the exponent 2(utv) = 4mt2 (m=0,1,..., (NHT-1)/2), 
the summation constraint becomes eq. 3.12; and when the exponent 2(utv) = 
4mt4 (m=0,1,...,(NT-2)/2), the summation constraint becomes eq. 3.13 
Since the summation must be non-negative when evaluated at s=jw. 


Collecting terms and simplifying, eqs. 3.12 and 3.13 become 


> 
r 5 (ape el) Sa Cues 
TYAS ll lia 
5 ec : ) 20 aS 
rq Le a eel - oy: 


u,r,q,m=0 
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u=0,1,...,(NT-1)/2 u=0,1,...,N4T/2 

r=0,1,...52 r=0n).,.. . 92 

q=0,1,...,T q=0,1,...,T 

m=0,1,...,N+T-1/2 m=0,1,...,(NT-2)/2 
(3.14) eels) 


and 


y = 4mt4-2u-4 e = 2u-q 


N = Order of original impedance. T order of augmenting polynomial. 
The final value of the indices is chosen so that if the value is non- 
integer the rule is to round down to the next lower integer. 
In the summations above, given values for a and 8, 6's must be chosen 
SO that these inequalities are satisfied. 


Expanding the generalized constraints for a first-order augmenting 


polynomial (i.e., T = 1; r,q = 0,1) 


2 2> 
Ad, + B65, + 66, = 0 
: ; (3.16) 
> 
DS, + E5,5, + F6, = 0 
where 
=: Opertl-2u°2utl 7 Oy ort? -2u-2u 
B= Ogm-2ue2ut) ~ mt? -2u°2u-1 
c= Oy m-2u-2u ~ “Amtl-2u82u-1 
| (3,17) 
DO Opeth 2uP2u ~ Camt3-2u°2ut] 
EE Cea ones CHpeoeBoe.1 
Ly 


- Oy rt3-2ue2uel _ Cpt? -2u°2u 
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[> 
ll 


jw 
il 


|= 
tl 


and 


GB) - By - AB, 


G38, + 0,8, + O48) 


lo 


O89 


O,By + 28, + OB, 


Ko) 


— OB, - OB, 


— OB, - AHR 


— ee 


0° 6 


- %6Bo 


(3.20) 


(322049 


(3.20b) 


(3.20c) 


so for a first-order augmenting polynomial the constraints are 


4456 


“Ket 


Us 


6 


2> 


+c OF 


0 


k 


k 


il 
© 
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ciemerk = O9F....,(N-1)/2 for the a's and k = 0.1,...,N/2 for 
the c's; N is the order of the original driving point impedance. 
Performing the expansions above for higher-order augmenting poly- 
nomials provides the results displayed in Table II. 
The equations of Table II can be put into matrix format. Equation 


3.21 becomes 








80 | bo 0 
k 1 
a 7 & = (3222) 
86 A aa 50 
= = 0 
0 ay Oy 
In general it can be shown that 
8° Qt 6 > 0 
a (3.22a) 
6° RK § > 0 
2 sj;2c 
where superscript t indicates a transpose, k = 0,1,...,N/2 and 


j = T+l for the first-order example (T is the order of the augmenting 


polynomial). 
ay 0 cy 0 ay 0 ; 
0 Ch 0 “ay 0 Chod 
co ae: ee = gr (3.23) 
0 “a, 0 Choy 0 “a 9 
zea ee etka? : 
x cleats Re “K-2 
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k+1 k k 
0 ab 0 “cy 0 ait 
-a 0 c 0 -a 0 
k k k-l = a (3.23) 
0 “ey 0 aed 0 “Cr 
oe : eT Soto 
0 ayy 0 “Cry 0 a9 


Note that Qu and RS are dimensioned (T+l) x (T+l). 
To illustrate the process so far developed consider the dp impedance 
used for the example in chapter V. 


a + 1/2s + l 
s +st1]1 


Z(s) = 


which requires a third-order augmenting polynomial. 


Ay = 1 a, = 12 dy = i 
Bo = 1 By = | Bo = | 


T+1=4 


us 
ll 


k = 0 for a's 
k = 0,1 for c's 


Ag = 0B, - By - A pBy = -3/2 


ay = 0 
Ne A8 g 2s 
7 C148» ot 


Zz 





-3/2 0O 
0 0 1 
Q, a 
-l 0 
0 S72 
l 0 
0 0 -3/2 
Ry a 
B72) 0 
0 -l 
It remains to solve 
where 


B. LINEARIZATION AND SOLUTION OF CONSTRAINTS 


-1 0 0 0 -1 
0 3/2 1 0 i 0 
0 0 2 -l 0 -3/2 
0 0 0 0 0 

| — 
3/2 0 iL 0 0 
0 -1 1 0 0 0 
R, = 

1 0 a! 0 0 i. 

0 0 0 -l1 0 


for the § coefficients according to eq. 


sl ile be aly 


(32225 


In this section the non-linear inequalities developed in the pre- 


ceding section are linearized using a set of congruent transformations 


and solved. 


becomes apparent. 


Note that they are symmetric for all orders. 


Inspecting matrices Q and R, the pattern of their formation 


Evalu- 


ating Q° and Re for fixed j=T+l over the indices k=0,1,...,N/2 produces 


a set of constraint inequalities as follow 
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\~ 

Lo 

|= 
lV 


(3.24) 


Kon 

v2) 

|~ 
lV 


or | Z 


|j~ 
pn 
Oo 
-) 
|\o~ 
V 
© 


ct 
i Lu. 
Kon 
lV 


(3.25) 


t ,N/2 
=; 


WV 


§ & 


Inequalities 3.25 represent N conditions involving T+l unknowns in 
i 

5. 6 ) 
Where 


—_— 


j-l 


As a first step toward the simultaneous solution of the set of 


inequalities above, consider the transforms for eq. 3.24. 
k 
wth i 
J (3.26a) 


6 = Ne x 
-"=-j- 


fox 
Cn 


where My Ww) is the normalized modal matrix of Qe (Ri) i.e., the 


matrix of normalized eigenvectors of gy (Ry). These orthogonal 


(1) 


Since one of the 6, can be chosen arbitrarily, there are actually T 
unknowns in 6. 
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transforms convert the quadratic forms 6° Q. (cima Ge. RY 6) into a linear 


combination of the squares of the coordinates of the vector x with no 


cross products. Alternately the transform can be viewed as reducing 


k 
the ‘Hus MAA R) to diagonal form; 


Ines), 
8° "6 ) (3.27) 
under the transform 
6 =M. x (3427 a) 
—j- 
becomes 
t k k _k fe k.-l1 _k .k 2 Z fi 
x (M.) Q,.M, x=x (M, pO oO Be ae oe ee Se Se oe 
x 8,7; * x @)) 8, =; * 00 ell joa 


(3228) 
: eed ou Kot k,-1 
since a property of an orthogonal matrix is Ct, ) = ore . The 
1's in eq. 3.28 are the eigenvalues of Q, and My is the orthogonal 


matrix of normalized eigenvectors which are always available. 


Using the transforms above in eq. 3.25 yields 


taelkwt kek. > 
x : .M. x = 0 
x Qt)’ gi Me x 
(3229) 
t kK Geek Kk > 
x N, R, N = 
7 N) wae 
Substituting 
k ou t k _k 
a : M 
q-J Qt) om ~J 
(3 .29a) 
px = cn‘) RX nk 
rj = J 7J 
where the D's are diagonal matrices, 3.29 becomes 
x" ps x 2 0 
= qj * 
(3.30) 
a pk x 5 0 
= rj = 
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Writing the quadratic forms as the product of a row and column vector 


eq. 3.30 becomes 


pk (row) z z 0 
TJ _ 
(3.31) 
k > 
ae (row) z = 0 
where 
2 
*0 
x 2 
1 
Z = ‘ (3 .3la) 
ol 


and p< (row) ¢ pé (row)) indicate p< ( p¥) Written as a row. Tiere 
q-J rj Qq-jJ tJ 
are N constraints on T+l variables in z. 

The solution to this linearized set of inequalities can be quite 
elegantly found using Danzig's simplex neenea- modified to permit 
optimization according to the criteria that as many of the individual 
inequalities as possible set to zero. The simplex method initially 
determines if the set of inequalities is consistent. For the purpose 
of this thesis it is recommended that one start with T = l. If the 
inequalities are inconsistent simply increase the order of the aug- 
menting polynomial by one and proceed until a feasible set is generated. 
According to Gievinimnaiies Pesitnegeven sparevor Z(s) has no jw-axis zeros a 
finite feasible set can always be found. 

A rigorous proof of this assertion is not available in the literature. 
Next the inverse transformation of eq. 3.26a must be accomplished in such 
a way that the resultant 6's are equal (i.e., 6 


§ § =5 15 ees). 


ono 


Zo 





re -snould be noted’ that” since Me and N. are different, different 
5's can be expected for the same x's. In Eq. 3.27 postmultiplying the 
Me QW) by 'D (an), which are arbitrary diagonal matrices, doesn't 
alter the congruency of the transform (i.e., the quadratic form in 6 is 
still transformed to a linear combination of the squares of the co- 
ordinates of the vector x). Since M (we) is the normalized modal 
matrix, postmultiplying by an arbitrary diagonal matrix is equivalent 
to scaling the modal matrix. 

Therefore, inverse transforming eq. 3.26a in such a way that the 6's 
k 


: : ae ; k 
are equal is equivalent to requiring that diagonal matrices a and a 


be found such that 
x (32-32) 


where eq. 3.32 is developed by setting the two equations 3.26a, 
postmultiplied by the required diagonal matrices, equal to each other. 

For the purposes of this thesis choose ws = ie I and we = oe: 
where and i are scalars and I is the identity matrix. With this 


substitution, eq. 3.32 becomes 


i* of x = mS né x (3538) 
=; = -j- 
or 
.x = BBE: né x (3 .33a) 
=j- -j = 
where 
k,k lee k 
Si 


Rewriting 3.33a by premultiplying both sides by wey" yields: 


wa! Me x= qk x (user) 
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Expanding eq. 3.33b over the index k,(j = T+l) yields 


tial 2 
cn’) mu? x = a x (n°) M. x = oO x (n’t y/ x = 0m 3. 
a os | a ae =a) — J 7 a i 
lett 0 1,0 Le cee jee lt “gNL2 1,N/2 

N. M, x = x N. M. x = x afar N, M, x = x 

(N,) mgs, dee x (N;) ls la: ae (N.) My #7 4 x 
(3.34) 

qv/? tut i Thee a?) "wx = gies. Cee lee * i ee 


Equations 3.34 are the familiar eigenvalue problem in which solution 
constants oN (the eigenvalues) can always be found for ws Me non- 
singular. Since equations 3.34 can always be solved for the eigenvalues, 
eq. 3.26a can always be inverse transformed such that the 6's are equal; 
and, consequently, a vector 6 can always be found which satisfies each 
of the constraints in eq. 3.25 simultaneously once a consistent set has 
been generated and solved for the vector x. 

Fortunately it is not necessary to solve for § explicitly. It is 


sufficient to realize that evaluating eq. 3.30 at the determined solution 


vector x yields numerical values: 


(3/35) 


If a particular scalar ky is zero, then the corresponding 5°Q55 


(O°R.5) is also zero; and if a particular ki, is positive, (the only 


teks) are also 


other possible condition) then the corresponding 5°95 (6 
positive although of different magnitude. This can be demonstrated as 
follows: 


Let 
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which is 
6= 1° Mex | (3Ng5e) 


fox ‘ 


Substituting eq. 3.35a in 3.24 yields: 


ee ee oe ee 2 k,2. 2 Kee 72 
x M) Gd) Q Mo x= 1) xg +1, yx +... +1, x,y 
k,2 2 2 2 
= (PY [oxy + 1x Fee FE] (3.35b) 


Comparing eq. 3.35b with 3.28, it is evident that eq. 3.35b is eq. 3.28 
multiplied by a positive constant which establishes the assertion above. 
Recalling that the quadratic forms are actually the numerator co- 
efficients of the real part of the augmented dp impedance, by solving 

eq. 3.31 and evaluating 3.30 the corresponding inequalities 3.24 are 

revealed as either zero or positive. It will now be demonstrated that 
this information is sufficient to determine the optimal augmenting and 
augmented polynomial to insure that the real part of Zz_(s) is positive 


term by term. 


C. OPTIMAL AUGMENTING AND AUGMENTED POLYNOMIALS 

In this section knowledge of the form of the augmented driving-point 
impedance (i.e., which coefficients are positive and which zero) will 
lead to the development of the optimal augmenting polynomial and con- 
sequently the optimal augmented driving-point impedance in a straight- 
forward way. | 


Consider an impedance function 


(3.36) 
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and an augmenting polynomial 


P(s) = ms + Dae (35364) 
Then 
~ ae ie ae — 
2 P(s) m," +n," set 


is the augmented polynomial. 


reeme 6 - Un.''n 
Rem, =) ———= pO eo (3.38) 
: (m,")* = (ny)? 
Z 2 
where 
ee ' ’ 
m, mim ate nj, 
= ' i 
Mo m)™m,) Fnon, 
aa | (3.39) 
gael hge li 
n= ' i 
No mn, + nom, 
Substituting eq. 3.39 in 3.38 yields 
' ' (aoe i = 
(m, +n, ) (m, ny ) (mm, nin) | 
ReZ (jw) = S=jw 
2 Co UR: ee re ey Co ee 
1 1 1 1 Z 2 (3.40) 
ReZ_ (jw) = Re(P(jw)) ReZ(jw) (3.41) 
P(jw) 
Num(ReZ, (jw)) a ReP( jw) Num(ReZ (jo) ) (3.42) 
Den (ReZ_ (jw) ) Den (ReZ_, (jw) ) ‘ 
where Num (ReZ , (jw) ) and Den (ReZ_, (jw)) represent numerator and 
denominator terms respectively. Therefore, 
Num (ReZ_, (jw) ) = ReP(jw) Num(ReZ (jw) ) (3.43) 


From eq. 3.24 we know which coefficients of Num(ReZ, (jw)) are zero 
and which are positive, and we also know Num(ReZ(jw)), the original 


unaugmented dp impedance. A simple-minded but norethe less effective 
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way to find ReP(jw) and the unknown coefficients of Num (ReZ , (j)) is 
to proceed by long division. The scheme is best illustrated by a 
numerical example. 
: : 8 6 4 2 
Given Num(ReZ(jmw)) =w - 20 + 3h - 4m + 5, developing the form 
of Num(ReZ_ (jw)) using the Danzig method outlined above to solve eq. 
3.31 for z and substituting into eq. 3.30, the following information 


is obtained concerning Num(ReZ, (jw)) 


Num (ReZ, (jw)) = Loa HO on nn ec sorte 6° 


where c, > 0 andc, > QO. 


5 6 
Dividing Num(ReZ(jw)) into Num (ReZ_, (jw) ) 
il +2 +1 
eeerses4e5. «(| I +0 +0 +0 +0 tee = 
1 -2 +3 -4 +5 
+2 -3 +4 -5 aheic 
+2 -4 +6 -8 +10 
rah] = +3 #(e,-10) te 
+1 -2 +3 -4 +5 
(c,-6) (c,-5) 


and setting ce = 6 and Si 5 the division terminates. The numerator 
of the ReP(jw) = a + NG + 1 and the numerator of ReZ , (jw) = us eS 
A feature of this method that can be useful is that the augmenting poly- 
nomial ReP(jw) always has non-negative coefficients. It may even be 
Hurwitz. In the Miyata method this added constraint is not needed; but 
in transfer function synthesis, where augmentation may be called for, 
the augmenting polynomial is additionally required to be Hurwitz. 

It turns out that if Z(s) is of second order or less, it is unnecessary 
to have exact knowledge of the form Num(ReZ_ (jw) ) since Poincare and Tegel 
have demonstrated that for Z(s) quadratic an augmented polynomial having 
but three non-zero positive terms can invariably be generated. For Z(s) 


of greater than second order, in general, exact knowledge of the form of 
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Num (ReZ, (jw) ) is required and the preliminary steps which yield this 
form cannot be bypassed. 

The problem is now essentially solved, and a brief outline of the- 
computational steps is in order. 

First generate a consistent set of non-linear inequalities (3.24) 
and then employ the transform 3.26a to set up and solve the resultant 
linear set of inequalities 3.31 by the Danzig simplex method which 
insures that as many as possible of the inequalities are zero, thereby 
minimizing the number of elements in the subsequent network synthesis. 
Armed with the solution vector x, determine the form of the solved 
inequalities 3.35 that yields which coefficients of the numerator of 
the real part of the augmented dp impedance are zero and which are 
positive. Divide the given Num(ReZ(jw)) into Num(ReZ_ (jw)) and solve 
for the unknown positive coefficients. The process is now complete and 
the numerator of the real part is known (evaluating the denominator is 
straightforward), and one can proceed to the synthesis part of the 
problem with a polynomial that is positive term by term. 

At this point it seems wise to demonstrate the solution algorithm 
and digital computer program developed to augment a moderately involved 


driving-point impedance. 
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IV. DIGITAL COMPUTER PROGRAM 


From consideration of the computational scheme outlined, it is 
obvious that the technique is admirably suited to implementation on a 
digital computer. This is espacially true since solving a set of linear 
inequalities under the constraint that the solution variables be greater 
than zero is very crater to the basic linear programming problem which 
is the subject of the Mathematical Programming Sere eons (MPS/360). 

Therefore, a highly developed and efficient system is readily 
available to solve the linearized inequalities, and the programming 
problem becomes one of simply generating the inequalities from the 
original driving-point impedance and inputing the set into the MPS 
system. 

The digital program required for systematic solution is straight- 
forward and, aside from typical programming difficulties, simple. The 
entire solution process has not been consecutively programmed due to 


interface restrictions between 0S/360 and MPS/360. The program follows 


in block diagram and Fortran format. 
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SUBROUTINE GROAN Calculate the element values for the 


=] 
x Z matrix of nonlinear inequalities. 
V 
SUBROUTINE GRUNT Form the matrix of inequalities. 
V 
SUBROUTINE GROWL A Convert the matrix of nonlinear 
inequalities to a matrix of linear 
inequalities. 
V 
SUBROUTINE PUNCH Prepare an input data deck for the MPS/360 
program which solves the linear inequalities 
and consequently the form of the original 
nonlinear inequalities. 
n=ntl 
V > 
MPS/360 Output the required order of augmentation 


and form of the numerator of the real part 


of the driving-point impedance. 
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V. SAMPLE PROBLEM AND SOLUTION 


5c + 1/2s + 1 


af + s+ l 


Z(s) = 


From computer program results 


NRez , (jw) = Tp eye Aan + Ose +c" te, 
: 4 2 
NRez(jw) =w - 3/2nH + 1 
DReZ(jw) = i ee + 1] 
Dividing NReZ , (jw) by NReZ(jw) gives: 
1 +3/2 +5/4 +3/8 
1 -3/2 + 1 1 0 0 0 te), tc, 
1 372, 1 
3/2 -1l 
3/2 -9/4 +3/2 
+5/4 -3/2 
+5/4 -15/8 +5/4 
+3/8 -5/4 
+3/8  -9/16  +3/8 
-11/16 -3/8 


Choosing Cy and c, equal to 11/16 and 3/8 respectively, the division 


3) 
terminates and the real-part numerator of the augmenting polynomial is 


seen to be 


Wore / Diets i aaa /8. 


The real-part numerator of the augmented polynomial has three non- 
zero terms and reads 
oy O411/ 160"+3/8. 
The denominator is 
DReZ (jw) x ReP(jw) 
which equals 


0 1 2843 (nS / Bw +7 (8243/8. 
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Therefore 


ie tena ye 


ReZ (jw) = oS 
Mee ih Pay cen 5 ey es -83/ 8 


which is positive for all w, as required. 
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TABLE ie 


COMPARISON OF THE NUMBER OF CIRCUIT ELEMENTS REQUIRED FOR SINGLE-N 
AND 1/2-N SPLIT MIYATA AND MODIFIED BOTT-DUFFIN METHODS 


Sing le-n Half-n Modified 
Split Spee Bott-Duffin 

eB wle a “le = “ie 23 
1 1 2 i 2 1 2 
2 4 3 4 3 S, 3 
3 9 4 6 4 7 S) 
4 16 5 10 ) 15 ii 
) 25 6 14 6 19 11 
6 36 7 17 7 33 15 
8 64 9 25 2 75 Syl 
10 100 11 i, 11 S13) 63 
12 144 13 44 13 915 127 
14 196 es) 52 15 635 255 
16 256 17 62 ly 1275 Seg! 
18 324 19 74 19 2550 1023 
20 400 21 86 21 5115 2047 
30 900 eB! n9 ou 163 ,835 65.955 


~ 
tt 


_ 
tt 


number of reactive components 


number of resistive components 
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